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Wilson's Numerical Renormalization Group method is 
used to study the paramagnetic ground state of the periodic 
Anderson model within the dynamical mean-field approach. 
For the particle-hole symmetric model, which is a Kondo in- 
sulator, we find that the lattice Kondo scale To is strongly 
enhanced over the impurity scale Tk\ Tq/Tk oc exp{l/3/}, 
where I is the Schrieffer- Wolff exchange coupling. In the 
metallic regime, where the conduction band filling is reduced 
from one, we find characteristic signatures of Nozieres ex- 
haustion scenario, including a strongly reduced lattice Kondo 
scale, a significant suppression of the states available to screen 
the f-electron moment, and a Kondo resonance with a strongly 
enhanced height. However, in contrast to the quantitative 
predictions of Nozieres, we find that the To oc Tk with a co- 
efficient which depends strongly on conduction band filling. 
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I. INTRODUCTION 

Lanthanide or Actinide-based heavy fermion com- 
pounds |j| can be viewed as a paradigm for strong cor- 
relation effects in solids. These systems comprise a 
great variety of different low-temperature states, includ- 
ing paramagnetic metals with either a Fermi liquid or 
non- Fermi liquid properties long range magnetic 

order, superconductivity or coexistent magnetism and 
superconductivity Q. While it is very interesting and 
challenging to understand the physics of these ordered 
ground states, not even the basic ingredient, namely the 
physics of the paramagnetic metal, has been fully cap- 
tured yet on a microscopic basis. The common under- 
standing is that the physical properties of heavy fermion 
compounds are largely determined by spin-flip scattering 
between spins localized on the Lanthanide or Actinidc 
sites and delocalized conduction quasi particles usually 
formed by the ci-states in the system. In the case of 
dilute compounds this spin-flip scattering leads to the 
well-known Kondo effect, i.e. the screening of the local- 
ized spins by the conduction electrons and the formation 
of a Fermi liquid with a low-energy scale set by the so- 
called Kondo temperature, Tk depending exponentially 
and non-perturbatively on the system parameters |3| . In- 
terestingly, the physics of the metallic phase of concen- 
trated systems can in some cases, e.g. CeAl 3 , CeB 6 or 
CeCu6 jjj , be at least qualitatively understood by a pic- 
ture of independent, but coherent Kondo scatterers, with 
the low-energy scale set by the Kondo temperature of the 



dilute system. However, for UPt3, URu2Si2 or Yb4As3 
|ij or the compound LiV204 recently characterized as 
heavy fermion system j^] there seem to exist two distinct 
energy scales; one high temperature scale, T K , describing 
conventional incoherent Kondo scattering, and a much 
smaller scale, To, also called coherence scale in litera- 
ture, which marks the onset of Fermi liquid formation 
characterized by quasi particles with strongly enhanced 
effective masses. 

This apparent discrepancy in the experimental situ- 
ation has triggered numerous theoretical investigations, 
from which two dominant competing scenarios have 
emerged. 

One scenario, based on a diagrammatic perturbation 
theory succeeded in approximately mapping the 

concentrated system onto a set of independent, effec- 
tive impurities which at low temperature form a coher- 
ent Fermi liquid state due to the lattice periodicity. The 
theory further predicts the existence of a single energy 
scale which uniquely determines the high energy proper- 
ties like Kondo screening as well as the formation of the 
low-tempcrature Fermi liquid. Most important, apart 
from moderate renormalizations which usually lead to 
an enhancement, this lattice scale is given in terms of 
the energy scale of the corresponding diluted system JtJ 
and the one particle properties close to the Fermi energy 
are well described by a picture of hybridizing bands, lead- 
ing to a density of states with a (pseudo-) gap slightly 
above the Fermi energy ||. The large effective masses 
can be accounted for by the observation that the Fermi 
energy lies in the region with flat bands (7||j|. A fur- 
ther theoretical study, based on a variational treatment 
of the system with a Gutzwiller wave function and em- 
ploying the Gutzwiller approximation, seems to support 
this result in the sense that the energy scale calculated 
is enhanced over the corresponding scale of the dilute 
system, i.e. there is quite likely only one relevant energy 
scale. However, this variational scale is more strongly 
renormalized |l0[ |. 

The other scenario is based on Nozieres argument 
that in concentrated systems there will generally not be 
enough conduction states available to screen all of the 
localized spins. This situation is engendered by the fact 
that only the band states within Tk of the Fermi surface 
can effectively participate in screening the local moments. 
The number of screening electrons can be estimated as 
Pq(0)Tk <§C 1, where /O§(0) is the conduction band density 
of states at the Fermi level. Thus, in a concentrated sys- 
tem one should encounter an intermediate temperature 
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regime where all band states available for screening are 
"exhausted" - from which the name "exhaustion physics" 
was coined - and the system starts to resemble an inco- 
herent metal where only part of the spins will be screened 
jn]]. However, since the Kondo screening clouds are not 
pinned to a particular site they can move through the 
system with an effective hopping matrix element and a 
residual strong local correlation, because there can never 
be two screening clouds on the same site. Based on these 
arguments Nozieres thus suggests, that at low temper- 
atures the system effectively behaves similar to a Hub- 
bard model with a small number of holes and strong local 
Coulomb repulsion. The remaining entropy can now be 
quenched by forming a state with either long-range an- 
tiferromagnetic order or short ranged antiferromagnetic 
fluctuations, which give rise to a small low-energy scale 
and a corresponding heavy Fermi liquid, too |l3f| . ft is im- 
portant to note that these phenomenological arguments 
by Nozieres do not make any reference to a particular 
model or specific parameter regime of a model like the 
periodic Anderson model; thus, if Nozieres arguments are 
correct, exhaustion physics should actually be the generic 
situation. 

Although both scenarios seem to be in accordance with 
some experimental facts, they apparently fail to capture 
the whole story. Moreover, in some cases they rely either 
on approximations that are difficult to control on a mi- 
croscopic basis or are entirely based on phenomenological 
arguments, as it is the case with Nozieres' point of view. 
In order to discriminate which ansatz is correct or outline 
the borders where each of the two scenarios may be valid, 
a more thorough microscopic study of the properties of 
the paramagnetic phase is highly desirable. 

Recently, Quantum Monte Carlo (QMC) simulations of 
the periodic Anderson model (PAM) within the Dynami- 
cal Mean-Field Theory (DMFT) pfQ have shown that 
when the conduction band is nearly half filled Jl5|] , there 
is a single scale consistent with the predictions of Rice 
and Ueda [flO|| , but obviously inconsistent with Nozieres' 
picture of exhaustion. On the other hand, when the con- 
duction band filling is significantly reduced, the states 
available for screening of the local moments appear to be 
depleted near the Fermi surface and the lattice coherence 
scale is strongly reduced from the corresponding impurity 
scale |T|,[l7]]. A protracted evolution of the photoemis- 
sion spectra |f8| , and transport |T^] are predicted and can 
be understood in terms of a crossover between the two 
scales. Although a quantitative relationship between the 
two scales could not be established, the reduction was as- 
cribed to exhaustion p8| . Nozieres subsequently argued 
that exhaustion should lead to a significant reduction of 
the lattice scale and predicted that To is at most Nd(0)T^ 

In this contribution we present results for the zero- 
temperature properties of the periodic Anderson model 
obtained within the DMFT in conjunction with Wil- 
son's numerical renormalization group (NRG) calcula- 
tions pi. We will show that for systems with approxi- 



mately half filled conduction band, one finds a single en- 
ergy scale, in accordance with the independent impurity 
picture and which qualitatively behaves as predicted by 
Rice & Ueda [[To); moreover, in this parameter regime we 
do not find any hint towards the occurrence of exhaus- 
tion physics; which raises the question to what extent 
Nozieres' phenomenological argument, that only a frac- 
tion pg(0)7k < 1 of the band states can contribute to 
the screening of the f-spins, is valid. 

In order to observe a behavior at all that is in at least 
qualitative agreement with Nozieres' exhaustion scenario, 
one actually has to strongly reduce the carrier concen- 
tration of the conduction band. For these low-carrier 
systems the physical behavior completely changes and 
can now be understood in the framework of exhaustion. 
However, a quantitative comparison of the lattice energy 
scale and the one of the corresponding impurity model is 
possible with our method, and indicates a strikingly dif- 
ferent relation between the two than recently predicted 
by Nozieres ]l9[ ] even in this parameter regime. 

The paper is organized as follows. In the next section 
we will introduce the periodic Anderson model and its 
treatment within the DMFT using NRG. In section [II 
we will discuss recent developments in understanding the 
physics of the PAM in some detail and provide a basis to 
interpret the NRG results presented in section [j^. The 
paper will be concluded by a summary and outlook in 
section fv]. 



II. MODEL AND METHOD 

A. Periodic Anderson Model 

The standard model used to describe the physics of 
heavy fermion systems is the periodic Anderson model 
(PAM) 



■£ V »<(cL/k«r + fc.' 



(1) 



In ([[]), cj^ creates a conduction quasi particle with mo- 
mentum k, spin a and dispersion £k, with ^2 £ k = £c 

k 

as its center of mass; creates an f-electron with mo- 
mentum k, spin a and energy e f and n- is the number 
operator for f-electrons at lattice site i with spin a. The 
localized f-states experience a Coulomb repulsion U when 
occupied by two electrons and the two subsystems are 
coupled via a hybridization Vk. Although it is in general 
a crude approximation we will assume for computational 
reasons Vk = const, for the rest of the paper. 

Given the usually rather complex crystal structure 
of heavy fermion compounds, the question to what ex- 
tent the simple model (Q), which describes orbitally 
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non-degenerate f-states, is appropriate. However, es- 
pecially in Ce-based systems the relevant configuration 
is f , whose multiplicity will in general be reduced to 
a Kramers doublet in the crystal field; the other crys- 
tal field multiplets are generally well separated from this 
ground state doublet The situation is more compli- 
cated in Uranium compounds, where with equal proba- 
bility a f 2 state can be the lowest configuration. In this 
case the assumption of a Kramers doublet is of course 
not justified. Nevertheless, even for these compounds 
the PAM seems to provide an at least qualitatively cor- 
rect description, so that we assume it to be the relevant 
model for those compounds, too. 

In contrast to e.g. the Hubbard model, where for one 
dimension an exact solution via Bethe ansatz is available 
and a combination of many-body techniques was suc- 
cessful to exploit the ground state properties exactly, no 
such exact benchmark is available for the PAM. The only 
limit where the physical properties are almost completely 
known is the impurity version of the model ([!]) (sin- 
gle impurity Anderson model, SIAM) [^) , where a single 
site with f-states exists which couples to the conduction 
states. In this limit one finds for |e/|/V, (e/ + U)/V ^> 1 
the mentioned Kondo behavior, namely a screening of 
the local moment by the conduction states with an en- 
ergy scale (Kondo temperature) 



exp 



/ = 2V 2 p c (0) 



Ef e f + U 



(2) 



with Po(0) the density of states at the Fermi energy of 
the conduction band and A a cut-off energy. Note, that 
7k depends exponentially on the system parameters, and 
especially non-analytically on V 2 . 

It is important to stress that the occurrence of the 
Kondo temperature is intimately connected to the ther- 
modynamic limit, where a (quasi-) continuum of band 
states exists at the Fermi energy. This property makes 
reliable calculations with techniques usually suitable to 
treat finite-sized two or three dimensional systems, like 
exact diagonalization or quantum Monte Carlo, of lim- 
ited value, since they can only handle small to moder- 
ately sized systems. For the quantum Monte Carlo an 
additional problem is the sign problem, that becomes 
increasingly serious as U and the system size increase. 
Direct perturbational approaches like e.g. FLEX are by 
construction restricted to small values of U and generally 
fail to capture even the basics of the Kondo physics. No- 
table exceptions occur again for the SIAM. Here, straight 
forward second order perturbation theory in U (not self- 
consistent) is able to reproduce at least qualitative fea- 
tures in the special case of full particle-hole symmetry 
PHp2 |. More recently, a novel ansatz using a perturba- 
tion theory about the unrestricted Hartee-Fock solution 
for the SIAM 123] has been shown to give even quanti- 



tatively correct results for spectral functions and energy 
scales of the SIAM in the Kondo limit 01. 



B. Dynamical Mean-Field Theory 

The other non-trivial limit, where an exact solution 
of the PAM (Q) becomes at least in principle possible is 
the limit of large dimensionality, where the dynamical 
mean- field theory (DMFT) becomes exact Jl|-0]. Here, 
the renormalizations due to the local Coulomb repulsion 
become purely local [Q, i.e. one obtains a one-particle 
self energy independent of momentum 



E / (k, z) -> S / (z) 



(3) 



This property may be used to remap the lattice problem 
onto an effective SIAM again 25 2q|. The non-trivial 
aspect of the theory comes about by the fact that the 
medium coupling to the effective impurity is a priori not 
known but has to be determined self-consistently in the 
course of the calculation P7|. In particular, the self- 
consistency relation for the PAM reads 



G loc (z) = I de 



1 

z-e f -K(z)-^(z) ' 



v 2 



(4) 
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iSIAM/ 



G mAM (z) 



with Po(e) the density of states of the bare conduction 
band in (|IJ). The second line in ([!]) defines a generalized 
hybridization function A(z), which implicitly depends on 

(z) and is thus modified from its non-interacting form 
by the presence of correlated f-electrons on other sites in 
an averaged way. 

The remaining task is to solve the effective SIAM de- 
fined by the set of parameters {sf, U} and the general- 
ized Anderson width 



r(w) = -$sm{A(uj + irj)} 



(5) 



Note that the self-consistency condition (^) requires the 
knowledge of a dynamical quantity, viz the one-particle 
Green function for all frequencies. This immediately 
rules out techniques like Bethe ansatz, since it is impossi- 
ble to calculate dynamical correlation functions with this 
method. To perform this task nevertheless, a variety of 
different methods have been developed and applied dur- 
ing the past decade: Quantum Monte Carlo simulations 
exact diagonalization (m), extended second order 
per turbation theory (iterated perturbation theory, IPT) 
[p^ , p9| , resolvent perturbation techniques 
moment approach and Wilson's NRG J31- 

In this contribution we used the last method, for the 
following reasons. First, it is tailored to capture the low- 
energy physics of the Kondo problem with high accu- 
racy. Second, it is able to identify exponentially small 
energy scales. Third, it is non-perturbatively and thus 




local- 
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also produces the correct dependence of the Kondo tem- 
perature on the parameters (see (§)) and, last but not 
least, it has only little numerical restrictions on the choice 
of model parameters. Together with past years develop- 
ments [3^j3^ | we are able to study the low-temperature 
properties of the PAM with the NRG with good accuracy. 



C. Details of Wilson's NRG 

Quite generally, the NRG is based on a logarithmic 
discretization of the energy axis, i.e. one introduces a 
parameter A > 1 and divides the energy axis into in- 
tervals [A-(" +1 ),A-"] for n = 0,1,..., oo |@. With 
some further manipulations fl33"f one can map the origi- 
nal model onto a semi-infinite chain, which can be solved 
iteratively by starting from the impurity and successively 
adding chain sites. Since the coupling between two adja- 
cent sites n and n + 1 vanishes like A~"/ 2 for large n, the 
low-energy states of the chain with n + 1 sites are deter- 
mined by a comparatively small number A sta tcs of states 
close to the ground state of the n-site system. In prac- 
tice, one retains only those Abates from the n-site chain 
to set up the Hilbert space for n + 1 sites and thus pre- 
vents the usual exponential growth of the Hilbert space 
as n increases. Eventually, after tt-nrg sites have been 
included in the calculation, adding another site will not 
change the spectrum significantly and one terminates the 
calculation. 

It is obvious, that for any A > 1 the NRG consti- 
tutes an approximation to the system with a continuum 
of band states but becomes exact in the limit A — > 1. 
Performing this limit is, of course, not possible as one 
has to simultaneously increase the number of retained 
states to infinity. One can, however, study the A- and 
Agates-dependence of the NRG results and perform the 
limit A — ► 1, Abates — + oo by extrapolating these data. 
Surprisingly one finds that the dependence of the NRG 
results on A as well as on the cut-off A sta tcs is ex- 
tremely mild, in most cases is a choice of A = 2 and 
7V statcs = 300 . . . 500 sufficient @. 

While the knowledge of the states is sufficient to cal- 
culate thermodynamic properties, the self-consistency (Q) 
requires the knowledge of the one-particle Green function 
or, equivalently the knowledge of the single-particle self 
energy E* (z). Since the NRG scheme works with a dis- 
cretization of the energy axis, it corresponds to a discrete 
system and by construction the Green function consists 
of a set of poles and an appropriate coarse-graining pro- 
cedure has to be applied. During the last 15 years con- 
siderable progress has been made to extract dynamical 
properties with the NRG, too, and it has been shown to 
give very accurate results also for e.g. dynamical one- and 
two-particle and also transport properties J35||56|]. For 
dynamical properties the NRG works best at T = 0, and 
various dynamic correlation functions can be calculated 
with an accuracy of a few percent. Although less well de- 



fined for finite temperatures, its extension to T > also 
shows very good agreement with exact results J32],[56| . It 
is quite remarkable as no sum-rules (Friedel sum rule, to- 
tal spectral weight etc.) must be used as input for these 
calculations. On the contrary, they can serve as an inde- 
pendent a posteriori check on the quality of the results. 
The first application of the NRG to the DMFT known to 
us is the work of Sakai et al. [ pTp2| where the symmetric 
Hubbard and periodic Anderson model in the metallic 
regime have been studied. In their papers, the authors 
point out some difficulties in the progress of iterating the 
NRG results with the DMFT equations, which are in our 
opinion largely related to the necessary broadening of the 
NRG spectra. 

As has been shown in the work of Bulla ct al. J|3| , these 
difficulties can be circumvented if, instead of calculating 
Gf(z) and extracting T,'(z) from it, one calculates the 
self energy directly via the relation 



(6) 



(7) 



Gt{z) 

Ft{z) = ((Uflfe,fl))» , 
which originates from the equation of motion 

zGl(z) = l + {{[f a ,H S j AM ],fi))(z) 

= l + (e f + A(z))Gl(z) + UFf( z ) 

with 

W N z - e k 

k 



Both correlation functions G£(z) and F£(z) appearing 
in (|J) can be calculated with the NRG and it turns out 
that the quotient of them gives a much better result for 
S^(z) than the use of G f (z) alone [33| . 

Let us note one particular problem in dealing with 
the PAM in the DMFT+NRG. As we will see in the 



beginning of section IV, the effective hybridization for 



the PAM in the particle-hole symmetric limit exhibits a 
pole right at the Fermi level (see Fig. || on page ||). It is 
clear that such a pole will lead to numerical difficulties. 
However, the NRG allows to deal with such a structure 
in a very efficient way, namely by including this single 
state represented by the pole into the definition of the 
"impurity" defining the beginning of the NRG chain. 



III. RECENT RESULTS 

Early studies of the PAM using the DMFT concen- 
trated on the particle- hole symmetric case e c = 0, 
2ef + U = 0, i.e. nf — 1 and n c = 1 |15 33j3^] . To solve 
the SIAM, the authors employed QMC §6|,|§ and NRG 
p2[ . Although for this particular situation the system is 



4 



a band insulator for symmetry reasons, one can extract 
a Wilson Kondo scale from the SIAM impurity suscep- 
tibility. Interestingly, this energy scale is enhanced with 
respect to the corresponding SIAM Kondo temperature 
in accordance with the results by Rice & Ueda JlO| , p^ |. 
In addition, at low temperatures the system can order 
antiferromagnetically |lq] . To avoid confusion and make 
the discussion more transparent we will use Tq and Tk 
from now on the distinguish the relevant lattice energy 
scale and the Kondo scale, respectively. 

One possibility to break the particle-hole symmetry is 
by depleting the band filling via changing e c , but keeping 
Uf sa 1, which has been done with QMC for U/V 2 » 4 
H[. The resulting phase diagram turns out to be 
quite interesting for several reasons. First, one finds a 
suppression of the antiferromagnetic order with decreas- 
ing n c and around n c — 0.5 a region with ferromagnetic 
order emerges | fL6| , which for n c < 0.5 can be attributed 
to a ferromagnetic effective RKKY-type exchange. Sec- 
ond, from studies of the resistivity and optical conduc- 
tivity jl7 18 1, one can infer that for n c < 0.8 the physics 
of the paramagnetic metallic phase drastically change. 
While for n c > 0.8 there seem to exist only one en- 
ergy scale Tq > Tk, the resistivity data suggest, that 
for n c < 0.8 the onset of coherence is marked by a tem- 
perature scale T <C Tk [|HJl^] . In addition, the spectral 
function drastically changes as one decreases n c and the 
effective Anderson width (||) , which for n c ps 1 has a peak 
at the Fermi energy, starts to develop a dip for to = 0. 
These observations were taken as fingerprints of Nozieres 
exhaustion scenario |[7|,[l8) ■ Further evidence for this in- 
terpretation comes from the work by Vidhyadhiraja et al. 
p8[ , which is based on 2 nd order perturbation theory in 
U (IPT). These authors, too, find a similar behaviour as 
function of n c . In addition, since their method allows to 
study T = 0, they could extract the energy scale T from 
their data. Interestingly, they found a relation between 
Tq and Tk of the form To oc (Tk)" for their results, which 
is precisley the behaviour predicted by the phenomeno- 
logical theory of Nozieres ITS]. However, we note that 
this particular relationship To cx T K in J|^] was found 
only if the width of the Kondo peak was used to esti- 
mate Tq. When the ratio of the masses was used instead, 
the relationship rather follows To oc Tk ■ 

The estimate of To coming from Gutzwiller ansatz cal- 
culations on the other hand apparently fail to capture 
the essential physics in this parameter regime, since they 
predict a ratio To /Tk > 1 for all values of n c |n]] . 

Although the results from QMC and IPT are at a 
first glance very convincing, there remain some prob- 
lems. First, both series of calculations were done with 
a comparatively small value of U/V 2 ~ 4 and a sys- 
tematic study of the behaviour as function of U /V 2 es- 
pecially for larger values is clearly necessary. Second, 
the results were either obtained by QMC 0,|l|l or IPT 
p8| . However, for large U/V 2 , the identification of ex- 
ponentially small energy scales with QMC is problematic 
due to its restriction to finite temperatures. The IPT on 



the other hand leads to ambiguous results as mentioned 
before; in addition, as a perturbational approach in U, 
it certainly cannot produce exponentially small energy 
scales. Thus, for a quantitative description of the low- 
temperature phase and especially a reliable calculation 
of the low-energy scale To, a non-perturbative technique 
at T = is necessary. 



IV. NRG RESULTS 

Such a method has become available recently by the 
application of Wilson's NRG §| to the DMFT 
which we use to study the paramagnetic phase of the 
PAM within the DMFT at T = 0. In order to per- 
form the energy integral in ([!]) analytically and to be 
able to make contact to earlier results we study a sim- 
ple hypercubic lattice in the limit of dimensionality 
d — ► oo. With the proper scaling jl3| this leads to 
pg(e) = exp (— (e 2 /t*) 2 ) /(t*^/¥) for the noninteracting 
DOS of the band states. In the following we use t* = 1 
as our energy unit. 

Let us start by briefly discussing the particle-hole sym- 
metric situation, i.e. 2e f + U = and e c = 0. Here, it 
is expected from symmetry and shown by calculations 
|]l5| , ^2 137 1 , that the concentrated system exhibits a hy- 
bridization gap. This is shown in Fig. [I], where the /- 
DOS A f {uj) = -±Sm{Gf(u + i5)} is plotted for the 
SIAM (dashed line) and the PAM (full line). The model 
parameters are U = 2 (i.e. £f — —1) and V 2 — 0.2. The 
result for the SIAM nicely shows the well-known struc- 
tures, namely the charge excitation peaks at w ~ ±(7/2 
and the Abrikosov-Suhl resonance (ASR) at the Fermi 
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FIG. 1. /-density of states Af(u) for SIAM (dashed line) 
and PAM (full line) in the particle-hole symmetric case 
2ef + U — and e c = 0. The model parameters are U = 2 
and V 2 = 0.2. The inset shows a blowup of the region around 
the Fermi level. 

level. This latter structure, which can be regarded as an 
effective local level, basically leads to the hybridization 
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gap in the PAM. The enlarged view of the region around 
the Fermi energy also shows, that the width of this hy- 
bridization gap and the width of the ASR are of similar 
order of magnitude. As we will dicuss in a moment the 
corresponding energy scale in the lattice is actually en- 
hanced over Tk, which sets the width of the ASR. This 
result can readily be anticipated from the self-energy of 
the /-states in Fig. |[ In the SIAM (dashed line) and the 
PAM (full line) one observes a nice parabolic extremum 



w 



-0.40 0.00 0.40 




Note that the above result, viz that the system is a Fermi 
liquid, is not a priori apparent from Fig. |l|, since the spec- 
trum represents an insulator. From the self-energy it is 
however evident that this insulating state can be triv- 
ially accounted for by using the picture of hybridizing 
quasi-particle bands, one of which is located at the Fermi 
energy and describes an effective f-level. The characteri- 
zation of the particle-hole symmetric system as band in- 
sulator is a priori not the only possible scenario. Indeed, 
for small V 2 and U —> oo an alternative route to an insu- 
lating state is via a Mott-Hubbard transition as found by 
Held ct al. for non-constant Vk |39|]. We did not observe 
such a transition for the case Vk =const. studied here for 
values of U as large as U = 10, but also cannot exclude 
such a scenario on the basis of the data available. 

It is also quite illustrating to have a look at the effec- 
tive Anderson width as defined by (||). This function is 
rather featureless for the SIAM (dashed line in Fig. ||), 
but exhibits a very pronounced structure near the Fermi 
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FIG. 3. Anderson width Y(w) for SIAM (dashed line) and 
PAM (full line) in the particle-hole symmetric case and pa- 
rameters according to Fig. The inset shows a blowup of 
the region around the Fermi level. 



FIG. 2. E-f {cj + iS) SIAM (dashed line) and PAM (full line) 
in the particle-hole symmetric case and parameters according 
to Fig. |lj From the real part (b) the Hartree term U/2 has 
been subtracted. The insets show a blowup of the region 
around the Fermi level. 

in Sym£^ (u>+i6) at the Fermi energy (see Fig. |]a), which 
is accompanied by a linear region in 5fteS^(cj + iS) in 
Fig. |^b; the slope is negative in both cases and definitely 
smaller for the PAM, as can be seen from the inset to Fig. 
||b. This means that, at least as far as the self energy is 
concerned, the system can be viewed as a Fermi liquid 
with a quasi-particle weight 



Z~ x = \- 



dbJ 



m* /to , 



(8) 



u=0 



where we introduced the notion of the effective mass 



level for the PAM (full line in Fig. g), namely a pseudo 
gap plus a strong peak right at to — 0; one can in fact 
show that the latter is a (5-peak. This <5-peak can be un- 
derstood as emerging from the sharp quasi-particle band 
with dominantly f-character at the Fermi level. 

The results by QMC and IPT suggest that one can ex- 
pect drastic changes in the physics of the model if one 
breaks the particle-hole symmetry. There are actually 
two possible routes to accomplish this goal; the first is 
to keep e c = 0, and thus n c w 1, but increase U so that 
2ef + U > 0. An example for the spectrum and the cor- 
responding hybridization function for £/ = — 1, U = 6 
and V 2 — 0.2 is shown in Fig. |[ where we plot the DOS 
Af(uo) (Fig. ||a) and the hybridization function T(uj) (Fig. 
^b). For these particular parameters the occupancies are 
n c 1 and n/ ~ 0.92. As in the particle- hole symmet- 
ric case one finds the characteristic three-peak structure, 



G 



again with a hybridizationgap in the DOS of the lattice 
(full line, sec inset to Fig. |ja). Note that this gap now is 



3 








-4 -2 2 4 6 8 



to 



(b) 




-2-10 1 2 



CO 

FIG. 4. DOS A f (uj) (a) and Anderson width T(u>) (b) for 
SIAM (dashed line) and PAM (full line) in the asymmetric 
case. The parameters are e c = 0, e/ = — 1, U = 6 and 
V 2 = 0.2. The insets show a blowup of the region around the 
Fermi level. 

located above the Fermi energy and its width is visibly 
larger than the width of the ASR in the SIAM; pointing 
again towards an enhanced energy scale for the lattice, 
which is also confirmed by an inspection of the self en- 
ergy. A replica of this (pseudo-) gap is also visible in 
the effective hybridization function of the lattice in Fig. 

, which in addition shows a pronounced peak slightly 
above the Fermi energy. The origin of this peak is the 
same as for rif = 1; only it now has to be damped due 
to the finite lifetime of the quasi-particles for uj > 0. 
It is quite interesting to note that the value of r(0) is 
actually reduced from its value in the SIAM, but the av- 
erage of T(u>) over a region of the order To around the 
Fermi energy is enhanced. If one assumes that it will be 
such an averaged value that determines the low-energy 
scale, one can readily understand that To is enhanced 
over Tk by this simple picture. We should mention that 
our results are very similar to those obtained more than 
10 years ago with resolvent perturbation theories 



In fact, the physical situation studied then was basically 
the same, namely a particle-hole symmetric conduction 
band hybridizing with an asymmetric /-level, however 
with U — oo. The interpretation of the structures in 
r(w), which may be viewed as the effective DOS of the 
medium visible to the impurity, then leads to the picture 
that due to the coherent admixture of /-states to the sys- 
tem, there will be an effective enhancement of T(u>) close 
to the Fermi energy. In this sense one may identify the 
physics in this region of parameter space with the picture 
of coherent Kondo scatterers. 

The most interesting question is how the energy scales 
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FIG. 5. Inverse effective mass as function of 1/7 (see text) 
for e c = and e/ = — 1. The inset shows the ration 
™pam/™siam together with a fit mp AM /mg IAM = e" 1/3/ . 

of the dilute system and the lattice, Tk and To, are 
related in this parameter regime. Let us recall that 
from the Gutzwiller ansatz one obtains Q T K /T = 
™pam/ to siam = exp(-l/27), where I = 8p^(0)V 2 /U; 
this prediction was found to be consistent with recent 
DMFT QMC simulations fl|] where the Wilson Kondo 
scale was estimated form the excess impurity susceptibil- 
ity. 

At T = 0, the most efficient way to extract the low- 
energy scale is by calculating the effective mass accord- 
ing to (JsJ) . With m/m* oc To we are then able to dis- 
cuss the variation of Tk and To with U. The result is 
shown in Fig. [|, where m/m* is plotted versus 1/7 and 
I = 2p c {Q)V 2 U/ (\e f \(e f + U)) is the Schrieffer- Wolff ex- 
change coupling. As already predicted from the spectra 
and self energies, the general relation To > Tk holds. 
In addition, both Tk and Tq apparently behave expo- 
nentially as function in 1/7, but with different slopes; 
in qualitative agreement with the predictions the slope 
of T) is smaller. To quantify the relation between Tk 
and T) further we show in the inset to Fig. [s] the ratio 
Tk/Tq — "ipAM/ m siAM as function of 1/7. Again, an 
exponential behavior is observed; however, in contrast to 
the predicted factor 1/2 we find ln(Tk/T ) = —1/37 from 
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FIG. 6. /-density of states Af(u>) (a) and effective hy- 
bridization function T(u>) (b) for PAM (dashed line) and 
SIAM (full line). The parameters are e c = 0.5, 2e/ + U = 0, 
U = 2 and V 2 = 0.2. 

our data. For the time being we do not have an explana- 
tion for this discrepancy between the QMC and NRG re- 
sults; whether it might be related to different schemes to 
extract To - scaling behavior of the excess susceptibility 
as function of temperature in the QMC versus effective 
mass at T — in the NRG - has to be clarified. 

The second possibility to destroy particle- hole symme- 
try is to choose e c ^ 0. As a typical result for that pa- 
rameter regime we show in Fig. || the DOS Af(uj) (Fig. 
^a) and the corresponding effective hybridization func- 
tion T(u) (Fig. |b) for e c = 0.5, 2e f + U = 0, U = 2, 
and V 2 = 0.2. For the filling we obtain rif « 1 and 
n c ss 0.6 for both the PAM (dashed line) and the single 
impurity Anderson model (SIAM, full line). As usual, 
one sees the characteristic structures, namely the charge- 
excitation peaks at u> w ±17/ 2 and the Kondo resonance 
at the Fermi level. However, in contrast to the results 
with s c = 0, we do not find any hint of a hybridization 
gap in the lattice DOS; in fact, the DOS looks pretty 
much like that of a conventional SIAM. The major dif- 
ference to the DOS of the SIAM is an enhancement of 
the ASR and a reduction of its width |17) , as is apparent 



from the inset to Fig. ||a. Particularly interesting is the 
behavior of the effective hybridization function in Fig. 
^b. Similar to the case e c = 0, 2ef + U > it is con- 
siderably reduced in the region around the Fermi level; 
in contrast to the former case, however, the sharp quasi- 
particle contribution is missing and the average value of 
r(<x>) over the region of the order Tq around uj = is still 
reduced from the non-interacting value here. The deple- 
tion in T(oj) around the Fermi level has been coined as 
hall-mark of exhaustion physics in the PAM and related 
models ]l7| , |38[ , since according to Nozieres phenomeno- 
logical picture jl9| the effective density of medium states 
available at a given site should be reduced due to screen- 
ing at other sites. 

The fundamental difference in the physics between 
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FIG. 7. c-density of states p c (u) of the PAM for e s = -1, 
V 2 = 0.2 and U — 6, e c = (n c w 1, upper panel) and {7 = 2, 
e c = 0.5 (n c ~ 0.6, lower panel). The dashed lines depict the 
bare conduction DOS Pq(uj) for the corresponding e c . 

n c fts 1 and n c <C 1 also manifests itself in the behavior 
of the DOS of the conduction states. Typical results for 
this quantity are shown in Fig. |7] for ej = — 1, V 2 = 0.2 
and e c = 0, U — 6 (n c w 1, upper panel) and e c = 0.5, 
U = 2 (n c w 0.6, lower panel). For comparison, the bare 
band DOS for the corresponding value of e c is also in- 
cluded (dashed curves); the choice of different values of 
U for the cases e c = and e c — 0.5 is necessary to ensure 
that both systems are metallic. One observes fundamen- 
tal qualitative differences in the DOS, especially close to 
the Fermi surface, which in our opinion are related to 
the different physical properties of the different regimes 
and do not depend critically on a particular choice of 
the interaction strength. For n c ~ 1, the appearance of 
a gap slightly above the Fermi level in the conduction 
DOS again supports the notion of hybridizing bands in 
this region of parameter space. On the other hand, the 
conduction DOS for the case n c « 0.6 does not show the 
typical form of hybridized bands, but merely a pseudo- 
gap at the Fermi energy. The fact, that the spectrum 
only develops a pseudo-gap with finite DOS for u) = 0, 
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can again be understood as sign of exhaustion, since the 
formation of a full hybridization gap requires complete 
Kondo screening by the band states to occur for each 
f-site; while the formation of a pseudo-gap can be inter- 
preted that only part of the f-sites are screened by the 
conduction electrons. 

The reduced effective hybridization around the Fermi 
level observed in Fig. ^| gives also rise to a reduced low- 
energy scale, characterized by an effective mass m* jm ~ 
17 in the PAM, whereas the corresponding quantity for 
the SIAM is m* /m as 8. The behavior of the low-energy 
scale as function of n r for fixed U = 2, V 2 = 0.25 and 
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FIG. 8. m/m* vs. n c for U = 2, V 2 = 0.25 and n f « 1. 
The lines are guides to the eye. The inset shows the ra- 
tio n^siAM/ m PAM = To/Tk versus n c . The data are fit by 
Tq/Tk oc n c exp (c ■ n c ) (full line), where c w 5/2. 

71/ sa 1 is shown in Fig. || Note that for n c 1 the 
value Tq is again enhanced over the impurity scale with 
m(T K /T ) = -1/37. Below n c « 0.8 the ratio T /T K 
decreases below one and falls rapidly and monotonically 
with n c , being almost two orders of magnitude smaller 
for n c w 0.2. The ratio "^siAM/ m PAM = ^b/^K is shown 
in the inset to Fig. ||. This ratio falls more rapidly than 
Tk, i.e. Tq/T k is not constant. The ratio Tq/Tk can 
be fit to a form Tq/Tk oc n c ■ exp (c • n c ) with c w 5/2 
(Full curve in the inset to Fig. |^). It gives an excellent 
account of the data. This behaviour especially means 
that Tq ~ n c T K asn c ^0@. Note that the Gutzwiller 
result for n c <C 1 predicts both Tq and Tk to behave like 
n c e _c '" c , but (i) predicts Tq/Tk > 1 and (ii) clearly gives 
no proportionality to n c in the ratio Tq/Tk- 

Nozieres phenomenological arguments also lead to 
an estimate of Tq as function of n c , namely Tq oc 



(t*) 2 /pS(o) 

with IPT EST 



19] . This relation has recently been tested 
and, within this approach, found to be 



fulfilled at least for U/V 2 m 4 between 0.4 < n c < 0.8 



when To is estimated from the width of the Kondo reso- 
nance. We first note, that clearly our result in Fig. || is 
not compatible with this prediction. In order to clarify 
the relation between Tq and Tk we compare Tq and Tk 
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FIG. 9. m/m* as function of U/V 2 for n c ~ 0.6 (squares 
and circles) and n c w 0.31 (diamonds) with nf ~ 1. The 
circles (squares) denote varying U (V 2 ) for fixed V 2 = 0.2 
(U = 2). The diamonds show the behaviour of m/m* for a 
different filling n c w 0.3 and V 2 = 0.2 for varying U. 

as function of U/V 2 for fixed n c . The results for n c sw 0.6 
as function of U/V 2 for both varying U at constant 
V 2 = 0.2 (circles) and varying V 2 at constant U = 2 
(squares) is shown in Fig. ^| on a semi-logarithmic scale; 
to study the dependence on the filling n c , we included 
calculations for n c — 0.31 (diamonds, varying U only). 
Evidently Tk and To both follow an exponential law 
T , T K oc exp(-o/7), where I = 8p^(0)V 2 /U is the 
Schrieffer-Wollf exchange coupling for 2ef + U = 0. How- 
ever, the curves for the SIAM and PAM for fixed n c are 
parallel in the semi-logarithmic plot in Fig. |^, i.e. the 
coefficients of U/V 2 in the exponents of both quantities 
are identical. This of course means To oc Tk rather than 
To oc (Tk) 2 , as predicted by Nozieres. Thus, neither as 
function of n c nor as function of U /V 2 does the lattice 
scale Tq obey Nozieres prediction. 



V. SUMMARY AND CONCLUSION 

We have presented results for the PAM obtained within 
the DMFT at T = using Wilson's NRG approach. 
For the range of parameters studied here, the system 
can always be characterized as a Fermi liquid with a 
strongly enhanced effective mass; this lattice scale To is 
enhanced over a corresponding impurity Kondo scale Tk 
for the particle-hole symmetric conduction band in ac- 
cordance with perturbational results or those from 







Gutzwiller ansatz |10(]. Moreover, the picture of hybridiz- 
ing quasi particle bands leading to (pseudo-) gaps in the 
DOS was found to be valid here. On the other hand, 
in the case of an asymmetric conduction band and espe- 
cially for low carrier concentration n c < 0.8, the spec- 
tral functions and corresponding effective hybridization 
functions show the signs of exhaustion and we observe 
a corrresponding strong reduction of the lattice scale To 

us 

Together with our results, an extensive picture of ex- 
haustion physics in the infinite-dimensional PAM has be- 
gun to emerge. Close to half filling the low-temperature 
properties of the model can be characterized by one en- 
ergy scale To > Tk; whereas away from half filling, two 
scales are apparent: Tk, where screening begins, and 
To <C Tk where coherence sets in. At low temperatures 
T sa To, the quasi particle features in the single-particle 
spectra become pronounced and have predominantly f 
character [Q. Since only the states near the Fermi en- 
ergy can participate in screening, the f-electron moments 
are screened predominantly by other f-electron states. At 
a temperature T m (T < T m < Tk), a dip begins to de- 
velop in the effective hybridization rate at the Fermi sur- 
face r(w w 0) fnl], indicating that the states available 
for screening are becoming exhausted. This is a direct 
confirmation of the qualitative features of Nozieres ex- 
haustion scenario. Nevertheless, for fixed n c , we find that 
To cx Tk, and as the conduction band filling changes we 
find T) oc n c exp(cn c ), both are in direct contradiction 
with the predictions of Nozieres. Thus, we conclude that 
Nozieres phenomenological arguments are too crude to 
capture the quantitative features of exhaustion. 

To our knowledge there has not yet been a systematic 
experimental study concentrating on the signatures of 
exhaustion in metallic HF compounds. However, there 
are several experimentally relevant consequences of ex- 
haustion predicted by DMF calculations. Most of these 
predictions are related to the presence of two relevant 
scales and the protracted behavior of measurements as 
a function of temperature in crossover regime between 
these scales JT(|. For example, as compared to the pre- 
dictions of the SIAM with the same To, the photoemis- 
sion peak should evolve much more slowly with temper- 
ature p8[ . In addition, it should have significantly more 
weight since the height of the peak goes like 1/T(0) and 
its width is set by Tq. Although these features have been 
reported inphotoemission experiments on Yb-based HF 
materials ^2) , these results remain controversial j43j and 
it has been suggested that the experimental spectrum is 
representative of the surface, and not the bulk |44| . For- 
tunately, transport and neutron scattering experiments 
probe much further into the bulk, and should display 
characteristic features due to exhaustion. The calcu- 
lated resistivity displays a non-universal peak and two 
other regions typical of HF systems and associated with 
impurity-like physics at high temperatures T rs Tk and 
Fermi liquid formation at low T < T [ pT[ . The peak 
resistivity shows features characteristic of exhaustion. It 



occurs at T w T m , the temperature at which exhaustion 
first becomes apparent as a dip T(tu sa 0). T m is non- 
universal in that it increases with decreasing n c , Tk and 
To. The Drude peak in the optical conductivity persists 
up to much higher temperatures than predicted by the 
SIAM and the Drude weight weight rises dramatically 
with temperature Jl7t . The quasielastic peak in the angle 
integrated dynamic spin susceptiblity also evolves more 
slowly with temperature than predicted by the SIAM, 
and it displays more pronounced charge-transfer features 

Despite the rich picture which has begun to emerge 
from DMF calculations, many questions remain unre- 
solved. Among these, three seem most prominent to us. 
First, it is unclear what is happening as n c — > 1. In 
this regime, following Nozieres argument, there should 
also remain too few states to screen the moments so the 
exhaustion scenerio would seem to apply, too; neverthe- 
less, the T) is enhanced relative to Tk and all effects of 
exhaustion vanish. It is tempting to attribute this van- 
ishing of exhaustion to another energy scale associated 
with the bare gap that appears in the spectra as n c — > 1. 
However, in recent calculations for a model with f-d hop- 
ping such that the hybridization Vk oc e^, where there is 
no conduction band gap when rif = n c = 1, in the regime 
of strong f-d hybridization Tq 3> Tk is again found, sug- 
gesting that there must be some more fundamental rea- 
son for the absence of exhaustion [j45f . Thus, these latter 
results together with the ones presented here, surely de- 
mand for a critical reinvestigation of the phenomenology 
of exhaustion. 

Second, thus far, all calculations are for the orbitally 
non-degenerate models. The effect of orbital degener- 
acy and crystal field splittings has yet to be determined. 
However, in the limit of infinite orbital degeneracy, the 
Kondo scale would seem to be unrenormalized. 

Third, following Doniach's arguments, RKKY interac- 
tion and superexchange will compete with Kondo screen- 
ing in the formation of the ground state. In the present 
work, we have explicitely concentrated on the paramag- 
netic state, i.e. these types of exchange do not enter the 
calculation. However, within the DMF it is possible to 
study the influence of RKKY or superexchange on the 
mean-field level by either looking at the susceptibility 
]l5| or allowing for a symmetry-broken state. Gener- 
ally, since the RKKY exchange grows like J 2 , and the 
Kondo scale is exponential in —1/J, Kondo screening 
and hence exhaustion is expected to be most pronounced 
when the Kondo exchange J ~ V 2 /U is large. How- 
ever, we have found that exhaustion can dramatically 
reduce the relevant low-energy scale; which may extend 
the region where the magnetic exchange dominates in the 
formation of the ground state. Thus a systematic study 
of the magnetic phase diagram as function of U/V 2 is 
clearly desirable. 
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